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Abstract In order to make asteroseismology a powerful tool 
to explore stellar interiors, different numerical codes should 
give the same oscillation frequencies for the same input phy- 
sics. Any differences found when comparing the numerical 
values of the eigenfrequencies will be an important piece of 
information regarding the numerical structure of the code. 
The ESTA group was created to analyze the non-physical 
sources of these differences. The work presented in this re- 
port is a part of Task 2 of the ESTA group. Basically the 
work is devoted to test, compare and, if needed, optimize 
the seismic codes used to calculate the eigenfrequencies to 
be finally compared with observations. The first step in this 
comparison is presented here. The oscillation codes of nine 
research groups in the field have been used in this study. 
The same physics has been imposed for all the codes in 
order to isolate the non-physical dependence of any possi- 
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ble difference. Two equilibrium models with different grids, 
2172 and 4042 mesh points, have been used, and the lat- 
ter model includes an explicit modelling of semiconvection 
just outside the convective core. Comparing the results for 
these two models illustrates the effect of the number of mesh 
points and their distribution in particularly critical parts of 
the model, such as the steep composition gradient outside 
the convective core. A comprehensive study of the frequency 
differences found for the different codes is given as well. 
These differences are mainly due to the use of different nu- 
merical integration schemes. The number of mesh points and 
their distribution are crucial for interpreting the results. The 
use of a second-order integration scheme plus a Richardson 
extrapolation provides similar results to a fourth-order in- 
tegration scheme. The proper numerical description of the 
Brunt- Vaisala frequency in the equilibrium model is also 
critical for some modes. This influence depends on the set 
of the eigenfunctions used for the solution of the differen- 
tial equations. An unexpected result of this study is the high 
sensitivity of the frequency differences to the inconsistent 
use of values of the gravitational constant (G) in the oscil- 
lation codes, within the range of the experimentally deter- 
mined ones, which differ from the value used to compute 
the equilibrium model. This effect can provide differences 
for a given equilibrium model substantially larger than those 
resulting from the use of different codes or numerical tech- 
niques; the actual differences between the values of G used 
by the different codes account for much of the frequency 
differences found here. 

Keywords Stars • Stellar oscillations • Numerical solution 
PACS 97.10.Sj • 97.10.Cv • 97.90.+j 



1 Introduction 

Asteroseismology is at present being developed as an effi- 
cient instrument in the study of stellar interiors and evolu- 
tion. Pulsational frequencies are the most important astero- 
seismic observational inputs. It is evident that a meaningful 
analysis of the observation, in terms of the basic physics of 
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stellar interiors which is the ultimate target of the investi- 
gation, requires reliable computation of oscillation frequen- 
cies for specified physics. This is a two-step process, involv- 
ing first the computation of stellar evolutionary models and 
secondly the computation of frequencies for the resulting 
models. Lebreton et al. (this volume) provide an overview 
of the tests of stellar model calculations. Here we consider 
the computation of the oscillation frequencies. 

An evident goal is that the computed frequencies, for 
a given model, should have errors well below the observa- 
tional error, which in the c ase of the CoRoT m ission is ex- 
pected to be belo w 0.1 uHz dBaglin et al .1120061) For th e Ke- 
pler mission (e.g., Christensen-Dalsgaard et al, 2007), with 
expected launch in early 2009, selected stars may be ob- 
served continuously for several years and errors as low as 
10~ 3 jllHz may be reachable, particularly for modes excited 
by the heat-engine mechanism. Since errors resulting from 
numerical problems are typically highly systematic, they 
may affect the asteroseismic inferences even if they are sub- 
stantially below the random errors in the observed frequen- 
cies. This must be kept in mind in the assessment of the es- 
timates of the numerical errors. 

During the last decades a lot of codes obtaining numeri- 
cal solutions of an adiabatic system of differential equations 
describing stellar oscillations have been developed. In order 
to ascertain whether any possible difference in the descrip- 
tion of the same observational data by different numerical 
codes is due to physical descriptions or to different numeri- 
cal integration schemes, the inter-comparison of these oscil- 
lation codes in a fixed and homogeneous framework is ab- 
solutely necessary. Some effort has been already done in the 
past but only regarding pairs of codes. Some codes have also 
developed a lot of internal precision tests. However, there is 
a lack of inter-comparison of a large enough set of codes. We 
aim in this study try to fix a set of minimum requirements for 
a code to be sure that any difference found is only due to a 
different physical assumption. 

Ideally, for a given model there should be a set of 
'true' frequencies with which the results of the different 
codes could be compared. This ideal situation could prob- 
ably be approximated by considering polytropic models for 
which it is relatively straightforward to calculate the equilib- 
rium structure with essentially arb i trary a ccuracy (see also 
IChristensen-Dalsgaard an d Mullari 1 1994) . In practice, the 
situation for realistic stellar models is more complex. Ow- 
ing to the complexity of the stellar evolution calculation the 
models are often available on a numerical mesh which is not 
obviously adequate for the pulsation calculation. The effect 
of this on the frequency computation depends on the detailed 
formulation of the equations in the pulsation codes. These 
formulations are equivalent only if the equilibrium model 
satisfies the 'dynamical' equations of stellar structure, i.e., 
the mass equation and the equation of hydrostatic support, 
and this is obviously not exactly true for a model computed 
on a finite (possibly even relatively sparse) mesh. One might 
define a consistent set of frequencies for a given model by 
interpolating it onto a very dense mesh and resetting it to 



ensure that the relevant equations of stellar structure are sat- 
isfied. The model is fully characterized by the dependence 
on distance r to the centre of density p and the adiabatic 
exponent F\ = (<9lnP/dlnp) a( j, P being pressure and the 
derivative being at constant specific entropy. Thus one could 
interpolate p(r) and F\ (r) to a fine mesh, and recompute the 
mass distribution and pressure by integrating the mass equa- 
tion and the equation of hydrostatic equilibrium. Frequen- 
cies of this model should then be essentially independent of 
the formulation of the oscillation equation and would pro- 
vide a suitable reference with which to compare other fre- 
quency calculations. Such a test may be carried out at a later 
stage in the ESTA effort. 

In Task 2, for now, we aim at testing, comparing and op- 
timizing the seismic codes used to calculate the oscillations 
of existing models of different types of stars. In order to do 
so we consider steps in the comparison by addressing some 
of the most relevant items that must be compared regarding 
the seismic characterization of the models: 

- Step 1: comparison of the frequencies from different 
seismic codes for the same model. 

- Step 2: comparison of the frequencies from the same 
seismic code for different models of the same stellar case 
provided by different equilibrium codes. 

- Step 3: comparison of the frequencies for specific pul- 
sators. 

The work presented here is mostly focused on step 1. 
At this step three different equilibrium models have been 
used. Two of them have been computed using CESAM 
dMorel and Lebretonl 120071) . with 902 and 2172 mesh 
points, and a third one with 4042 mes h points provided 
by ASTEC dChristensen-Dals gaardl l2007l) .We present inter- 
comparisons using the two models with the larger numbers 
of mesh points. The same physics and physical constants 
(except G) are used for all the oscillation codes. Frequen- 
cies in the range of [20,2500] jxHz, belonging to spherical 
degrees £ = 0,1,2 and 3 have been calculated, in order to 
recover most of the possible values we can find in the obser- 
vational photometric data. 

We present in Sect. 2 the equilibrium models used and 
we analyze their main features. Section 3 is devoted to the 
different oscillation codes used, discussing common require- 
ments and variations in the treatment. In Sect. 4 the direct 
comparison of frequencies for different frequency ranges 
and spherical degrees is presented in detail. Sections 5 and 
6 analyze the same inter-comparison for the values of the 
large and small separations, respectively. Section 7 discusses 
the dominant effects that contribute to the differences in fre- 
quencies determined by the different codes. Conclusions are 
given in Sect. 8. 



2 The equilibrium models. 

To ensure that any difference obtained in the inter-compari- 
sons is only due to differences in the numerical schemes, 
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Table 1 General characteristics of the models used for the inter- 
comparison. 
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we imposed to all the codes the use of the same equilib- 
rium models. These models were supplied in several for- 
mats: DSC, FGDNG, SROX, and FAMDL. The first model was 
required to have 900 mesh points, and it was provided by 
CESAM. The differences in this case reached unacceptable 
values of 2-3 /iHz when the same integration schemes are 
compared, and even 10 /J Hz when the use or not of the 
Richardson extrapolation are compared. The maximum dif- 
ference found for large separations with this model is 1 /iHz. 
This showed that either a larger number of mesh points or 
Richardson extrapolation was necessaryQ 

The second model was also provided by CESAM (re- 
ferred from now on as M2k). It uses a grid, with 2172 
mesh points, more suitable for asteroseismic purposes. Gen- 
eral characteristics of the model are presented in Table 1. 
These are typical of a 7 Doradus star showing oscillations 
in the asymptotic g-mode regime, and also around the fun- 
damental radial mode. A priori, solar-like pulsations cannot 
be excluded for this type of star and therefore it can be a 
good candidate for a global study. In Fig. Q] A* (which is a 
quantity directly related to the Brunt- Vaisala frequency N 2 : 
A* = rg~ l N 2 , where g is the gravitational acceleration) is 
depicted as a function of the relative radius (x = r/R) in a 
region of steeply varying hydrogen abundance, and hence 
mean molecular weight ju, just outside the convective core. 
The model is in a phase of a growing convective core. If 
diffusion and settling are neglected this leads to a discon- 
tinuity in the hydrogen abundance and hence, formally, to 
a delta function in A*; also, there is a region of 'semicon- 
vection' at the edge of the core. In fact, the figure shows 
that three points in this model display an erratic variation in 
A* just in the transition region between the convective and 
the radiative zone; also, the mesh resolution of this region 
of rapid variation seems inadequate. This is emphasized by 
Fig. |2] which shows the distribution of mesh points of this 
model along the stellar radius, indicating that there are not 
enough mesh points in this transition zone. As discussed be- 
low, these features of model M2k give rise to frequency dif- 
ferences in the comparison, particularly for those modes for 
which this inner part is critical for their physical description. 

The third and last model is a 4042 mesh points model 
(from now on M4k) provided by ASTEC. General character- 
istics of the model are presented in Table 1 . It has been com- 
puted to have overall characteristics similar to the previous 
CESAM model in order to understand better the differences. 
However, as discussed by Christensen-Dalsgaard (this vol- 

1 Detailed results of this investigation can be found at 
http : //www. astro .up.pt/corot/compf reqs/task2/ 
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Fig. 1 A* (related to the Brunt- Vaisala frequency) as a function of the 
relative radius for the two equilibrium models discussed in the text 
in the ,u -gradient zone, close to the convective core. The mesh points 
provided in the models are indicated. 
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2 Accumulated number of layers as a function of the relative ra- 
for the two equilibrium models explained in the text. 



ume) particular care has been taken in the treatment of the 
jj. -gradient region; the semiconvective region was replaced 
by a region with a steep gradient in the hydrogen abundance, 
defined such as to ensure neutral stability of the temperature 
gradient. As a result, A* for this model, also presented in 
Fig. [TJ shows a well-defined and reasonable behaviour. The 
distribution of the mesh points can be found in Fig. |2j In 
the central and outer parts of the stellar model, the distribu- 
tion is similar to the CESAM model M2k. It is in the inner 
zones, particularly in the boundary region between the con- 
vective core and the fi -gradient zone, that the models present 
different mesh-point distributions, with M4k providing a far 
superior resolution of this critical region. 
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3 Oscillation codes and requirements. 

All oscillation codes involved in this task were asked to pro- 
vide adiabatic frequencies in the range of [20, 2500] fx Hz and 
spherical degrees I = 0, 1,2 and 3. In addition, the solution 
of the equations must satisfy the following requirements: 

- To use the mesh provided by the equilibrium model, no 
re-meshing is allowed. 

- To set the Lagrangian perturbation to the pressure to zero 
(dP — 0) as outer mechanical boundary condition. 

- To use the physical constants prescribed in Task 1 . 

- To use linear adiabatic equations. 

Nevertheless, some other schemes for the numerical so- 
lution of the differential equations (from now on called for 
simplicity "degrees of freedom") remain open. Nine oscilla- 
tion codes of different research groups in the field have been 
used in this inter-comparison exercise. A summary of the 
participating codes and the different "degrees of freedom" 
provided by each one is found in Table 2 and include: 

- Set of eigenfunctions: Use of the Lagrangian or the Eu- 
lerian perturbation to the pressure (dP or P'). This ob- 
viously affects the form of the equations; in particular, 
when using dP the equations do not depend explicitly 
on A*. 

- Order of the integration scheme: Most of the codes use 
a second-order scheme, but some of them have imple- 
mented a fourth-order scheme. 

- Richardson extrapolation: Some of the codes using a se- 
cond-order scheme have the possibi l ity to use Richard- 
son extrapolation (Shiba hashi et all 1 198 lh to decrease 
the truncation error; combining a second-order scheme 
with Richardson extrapolation yields errors scaling as 

.jV being the number of mes h points (e.g. 
IChristensen-Dalsgaard a nd Mullanl [T994l) . 

- Integration variable: Two integration variables are used: 
1) the radius (r), or 2) the ratio rjP. The latter variable 
may reduce the effect of rounding errors in the outer lay- 
ers (see Sect. |7.3l l. 

- Despite the requirement that the physical constants be 
fixed at the values for Task 1, the codes used slightly dif- 
ferent values of the gravitational constant G, as listed in 
Table [2] Ideally the equilibrium model should have been 
computed with the prescribed Task 1 value (6.6716823 x 
10~ 8 cgs) which should then have been used for the os- 
cillation calculations. In practice Model M4k was com- 
puted with G = 6.67232 x 10~ 8 cgs. Using different val- 
ues of G in the oscillation equations clearly gives rise to 
inconsistencies, with potential effects on the frequencies, 
as discussed further in Sect. 17. 11 



Note that most of the oscillation codes put I = in the 
general non-radial differential equations to obtain the ra- 
dial modes, except for LOSC that uses the LAWE differ- 
ential equation (Linear Adiabatic Wave Equation for radial 
modes), and for GraCo for which results will be shown for 
both sets of equations. It should be noted that the LAWE 
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Fig. 3 Frequency comparison (reference line is GraCo) for modes with 
I = obtained for the model M2k. ADIPLS and NOSC frequencies 
have been obtained without using the Richardson extrapolation. 



does not depend on A*. All the main characteristics and nu- 
merical schemes are presented in previous chapters of this 
volume. 

The Nice code (NOSC) has the options of using P' or dP 
as dependent variables, and r or r/P as independent variable. 
However, all NOSC results presented here use P' and r. 

The use of the mesh provided in the equilibrium model, 
rather than meshes optimized for the different kinds of 
modes, may result in inadequate resolution of the rapid vari- 
ation of high-order p- and g-modes and hence larger trun- 
cation errors in the solution of the differential equations. In 
these cases, therefore, the second-order scheme may result 
in unacceptable errors. Here the use of higher-order inte- 
gration schemes (a fourth-order scheme or a second-order 
scheme followed by Richardson extrapolation) is therefore 
expected to give better results. For low- and intermediate- 
order modes we expect little effect of the use of higher-order 
schemes. 



4 Frequency inter-comparisons. 

In this section, the results of the direct frequency inter-com- 
parison are presented. We have structured this study split- 
ting the frequency range in three parts for the non-radial 
case, only one for the radial case, and comparing codes only 
with similar selections of "degrees of freedom". In addition, 
the influence of using different selections in the same code 
(GraCo in this case) is shown; for these tests the value of G 
in the oscillation calculations was the same as was used to 
compute the equilibrium model. 



4.1 Radial modes. 

These are shown in a single frequency range. In Fig. [3] 
the results obtained using the model M2k are presented. 
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Table 2 List of participating codes in this inter-comparison and "degrees of freedom" for each code: (a) Eigenfunctions, (b) Order of the 
Integration Scheme (2 or 4), (c) Use (y) or not (n) of the Richardson extrapolation, (d) integration variable used, and (e) the choice of the 
gravitational constant G. The references for each code are also given. 
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The reference line for all the inter-comparisons is selected 
to be GraCo. In this figure the reference frequencies have 
been obtained using P', second order, no Richardson ex- 
trapolation and r as independent variable (see Table 2). 
Two sets of codes can be identified in the figure, ADIPLS- 
NOSC-GraCo, with differences lower than 0.25 /iHz, and all 
with the same "degrees of freedom", and OSCROX-PULSE- 
LOSC-POSC with differences for high frequencies lower than 
2 /iHz with GraCo, but with differences among them around 
0.5 /iHz. This second set of codes differs from the first one 
in the use of a fourth-order numerical scheme instead of a 
second-order one. In the present figure, in Fig. [5] and in Figs 
[7] and [8] showing inter-comparisons with model M2k, POSC 
has been chosen as representative of the codes using sec- 
ond order plus Richardson extrapolation (see Table O. We 
can see how the use of this integration procedure provides 
similar results as the fourth-order solutions. These differ- 
ences, around 0.5 /iHz for codes using the same integration 
scheme, and 2 /iHz for codes using different schemes, are 
larger than the expected precision of the coming observa- 
tional data. Therefore this effect can change any detailed 
physical description as interpreted by different oscillation 
codes. Also, we point out that the differences between the 
codes using fourth-order schemes and GraCo results using a 
second-order scheme indicate that the model has an insuffi- 
cient number of mesh points for asteroseismic studies. 

Figure [4] shows the results obtained for the model M4k. 
In the top panel models with second order, no Richardson 
extrapolation, r, are depicted. All the differences are lower 
in magnitude than 0.014/iHz, i.e., two orders of magnitude 
lower than those obtained for model M2k. Therefore im- 
proving the mesh, including a doubling of the number of 
points, provides a very substantial improvement in the pre- 
cision, making these values more acceptable for theoretical 
modeling. 

The middle panel of Fig. [4] presents the differences ob- 
tained with models providing fourth-order integration so- 
lutions or second-order plus Richardson extrapolation. The 
global precision here is similar to the previous case or even 
slightly lower, with differences lower than 0.02 /iHz. It is in- 
teresting to point out that we cannot directly distinguish be- 
tween a fourth-order integration scheme solution (OSCROX) 



or a second-order plus Richardson extrapolation (the rest). 
However, it is noticeable that the GraCo-OSCROX-POSC fall 
in one group and ADIPLS-NOSC in a second, with a slight 
difference in the latter case. These two groups are distin- 
guished by the value of G (cf. Table |2), with ADIPLS and 
NOSC having similar but not identical values. This pattern 
will be found in other cases also. The LOSC behaviour is dis- 
cussed in the next paragraph. 

Finally in the bottom panel of Fig. [4] a comparison be- 
tween the different "degrees of freedom" using only the 
GraCo code are presented. As reference we have used the so- 
lution with "degrees of freedom": X =(i = 0, no Richardson 
extrapolation, P', r). For each comparison we have changed 
only one of these "degrees of freedom" at a time, keeping 
the rest unchanged (solutions X'). The most prominent ef- 
fect arises from the use of Richardson extrapolation, which 
changes the frequencies by nearly 0.8 /iHz for the highest- 
order modes, substantially more than the expected observa- 
tional accuracies. For model M2k (see Fig. [3), we have sim- 
ilarly found a change of 2 /iHz, reflecting the smaller num- 
ber of mesh points. This clearly shows that second-order 
schemes are inadequate, even for the mesh in M4k, for the 
computation of high-order acoustic modes; as expected the 
effect decreases rapidly with decreasing mode order. The use 
of r/P as integration variable provides small differences, al- 
ways lower than 0.008 /iHz. These differences are of the or- 
der of those obtained in the top and middle panels. The use 
of the Lagrangian perturbation to the pressure as variable 
(dP) or the use of the LAWE differential equation provide 
very similar differences, lower than 0.05 /iHz, but with an 
oscillatory pattern. This pattern is very similar to that ob- 
served for LOSC, which also uses LAWE to obtain the ra- 
dial modes. As discussed in Sect. 17.21 this oscillatory pat- 
tern arises from an inconsistency in the thermodynamics of 
model M4k which affects A*; solutions of equations that do 
not depend on A* (i.e., the LAWE or the equations based on 
dP) are insensitive to this effect. Therefore, even for model 
M4k and radial modes, the use of different integration pro- 
cedures can give different values for the oscillation eigenfre- 
quencies that are non-physical in nature. These non-physical 
sources of differences are mainly some inconsistencies in 
the equilibrium models (see Sect. 17.21) and the lack of mesh 
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Fig. 4 Frequency comparison (reference line is GraCo) for modes with 
I = as a function of the frequency obtained for the model M4k. In the 
top panel the models with: second order, no Richardson extrapolation, 
r, are depicted (NOSC-ADIPLS-GraCo-FILOU-POSC). The middle 
panel presents the differences obtained for models providing fourth- 
order integration solutions or second-order plus Richardson extrapola- 
tion (LOSC-OSCROX-NOSC-ADIPLS-GraCo-POSC). In the bottom 
panel a comparison of the different "degrees of freedom" only using 
GraCo is presented. 



points. But when the same numerical schemes are used, the 
different codes provide very similar frequencies. 

4.2 Non-radial modes with 1 = 2 

To illustrate the differences appearing in the case of non- 
radial modes, the spherical degree I = 2 has been chosen. We 
have divided the frequency spectrum into three regions: 1) 
Large frequencies ([500,2500] /iHz), 2) frequencies around 
the fundamental radial mode ([80,500] ;UHz), and 3) low- 
frequency region ([20,80] jiHz). In all cases a study similar 
to that developed in the radial case has been carried out. 

For the sake of simplicity the high-frequency differences 
are not represented since the results are very similar to those 
presented for the radial case. Only a slightly higher preci- 
sion is found in this case. The results of LOSC present the 
same pattern as the radial case, owing to the use of 8P as 
eigenfunction in that code (see Sect. I7.2I) . 

Figure [5] shows the results obtained for model M2k 
when comparing 1 = 2 frequencies around the fundamen- 
tal radial mode. The main differences are smaller than in 
the high-frequency region, corresponding to the low order 
of the modes and the consequent lesser sensitivity to the 
number of mesh points. The largest differences are found 
for two modes of frequency near 345 and 362 /iHz show- 
ing avoided crossing; these modes have a mixed charac- 
ter, with fairly substantial amplitude in the fi -gradient zone. 
PULSE and LOSC present differences around 3 - 4 JiHz, POSC 
around 2/iHz, OSCROX less than 1 /iHz, and the rest do not 
present significant differences for these two modes. That is, 
the largest differences are found in these codes when us- 
ing a fourth-order integration scheme or a second-order plus 
Richardson extrapolation. The values of these differences 
are larger than the expected precision of the satellite data 
to come. They clearly reflect the inadequate representation 
of the fi -gradient zone in M2k, with higher-order schemes 
being more sensitive to the resulting inconsistency in the 
model. 

In the top panel of Fig. [6] the inter-comparisons of the 
frequencies (£ = 2) for model M4k, using a second-order 
scheme without Richardson extrapolation, are depicted for 
the same frequency range as in Fig. [5] The main differ- 
ences are two orders of magnitude lower than those ob- 
tained for model M2k, and they are also of the same or- 
der of magnitude as those obtained for the high-frequency 
range with M4k. All these differences remain always lower 
than 0.025 jtiHz. GraCo and POSC are extremely close, while 
FILOU and ADIPLS provide very similar results, with slightly 
larger differences for NOSC, by about 0.005 /iHz, relative to 
these codes. Thus again the differences are related directly 
to the different values of G. There remain small wiggles for 
the mixed modes near 350/iHz but reduced more than two 
orders of magnitude relative to the largest differences for 
model M2k, reflecting the superior resolution of the critical 
region in model M4k. 

The middle panel of Fig. [6] shows the same inter-compa- 
rison as the top panel for codes using a fourth-order scheme 
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Fig. 5 Frequency comparison (reference line is GraCo) for modes with 
I = 2 for the model M2k around the fundamental radial mode. ADIPLS 
and NOSC frequencies have been obtained without using Richardson 
extrapolation. 



or a second-order plus Richardson extrapolation. The preci- 
sion is similar to the previous case, with differences of the 
same order of magnitude. We can distinguish two groups of 
codes providing very similar results: OSCROX-POSC-GraCo 
and LOSC-ADIPLS-NOSC. This distribution does not depend 
on the integration scheme, but again reflects the values of 
G. The wiggle of the mixed modes is also very similar to 
that obtained with a second-order scheme without Richard- 
son extrapolation. 

Finally, an inter-comparison of different "degrees of 
freedom" using only GraCo is presented in order to test the 
differences obtained for the different choices. This compari- 
son is depicted in the bottom panel of Fig. [6] The use of the 
Richardson extrapolation is not very important, as expected 
for these modes of low order, with effects generally smaller 
than 0.002 /iHz, although a larger value is present for the 
mixed modes. Using dP as variable gives differences larger 
than the differences among codes with the same "degrees of 
freedom". The use of the integration variable r/P does not 
introduce significant differences. 

Compared with model M2k, we find a general reduction 
of the differences for M4k, the main effects being in the re- 
gion of mixed modes with an improvement reaching up to 
three orders of magnitude. This is obviously not a simple 
consequence of the doubling of the number of mesh points. 
The main reason is likely the inadequate resolution shown 
by model M2k in the description of the Brunt- Vaisala fre- 
quency in the region close to the boundary of the convec- 
tive core, which is not present in model M4k. The avoided- 
crossing phenomenon and the behaviour of the mixed modes 
are very sensitive to the detailed treatment of this region, in- 
cluding the effects of semiconvection. Therefore, an accu- 
rate description of N 2 is critical for the oscillation codes in 
order correctly to obtain eigenfrequencies for modes near an 
avoided crossing. 

The direct frequency inter-comparison ends with the stu- 
dy of the low-frequency region; as above we concentrate on 





r 




-0.005 r 


N 
I 




3 


-0.01 - 


> 




O 




o 

CO 


-0.015 - 


a 




> 






-0.02 - 




-0.025 - 




0.005 r 




- 




-0.005 - 


I 


3 




> 


-0.01 - 


o 
O 




CD 




6 
^> 


-0.015 - 




-0.02 - 




-0.025 - 




0.008 r 




0.006 - 




0.004 - 


X 




3 






0.002 - 


> 










- 




-0.002 - 




-0.004 - 




100 150 200 250 300 350 
Frequency QiHz) 



400 450 500 



Fourth order, 


second order + 


Richardson extrapolation 








i • ••**• « ». 


-W ....» „. 


• : * /'..'•* 






'■V -■' 






V/ / 


ADIPLS 
NOSC 
POSC 
LOSC 
OSCROX 







100 150 200 250 300 350 400 450 500 
Frequency (uHz) 









Rich » 








SP —— 








r/P a 




























:.7.--.S-^:;:v:g' 1 ' 































100 150 200 250 300 350 
Frequency (p.Hz) 



400 450 500 



Fig. 6 Frequency inter-comparison (reference line is GraCo) for 
modes with I = 2 as a function of the frequency obtained for the 
model M4k and for modes around the fundamental radial. In top panel 
models with: second order, no Richardson extrapolation, and r are de- 
picted (NOSC-ADIPLS-GraCo-FILOU-POSC). Middle panel presents 
the differences obtained for the models providing fourth-order integra- 
tion solutions or second-order plus Richardson extrapolation (LOSC- 
OSCROX-NOSC-ADIPLS-GraCo-POSC). In bottom panel an inter- 
comparison for the "different degrees" of freedom only with GraCo is 
presented. 
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Fig. 7 Frequency comparison (reference line is GraCo) for modes 
with i = 2 obtained for the model M2k in the low frequency region. 
ADIPLS and NOSC frequencies have been obtained without using the 
Richardson extrapolation. 
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Fig. 8 Period separation as a function of the frequency in the asymp- 
totic g-mode region for model M2k. 



modes of degree 1 = 2. Fig. [7] shows the results obtained for 
model M2k. In this case the differences are lower than ob- 
tained at higher frequencies for this model (cf. Figs [3] and 
[5). The most surprising behaviour is present in PULSE and 
LOSC. POSC and OSCROX also present some differences in 
the same region. To understand the reason for these differ- 
ences and the region where they appear Fig. [8] shows the 
period separation ATI between adjacent modes. The first- 
order asymptotic g-mode theory predicts a constant separa- 
tion of the periods in this regime for a given I. However, 
when the equations are solved numerically, this period spac- 
ing presents several minima, a nd these minima are di rectly 
linked with the mode trapping dBrassard et all fl992bh . Fig. 
[8] shows that the position of one of these minima is the same 
as the position of the largest differences. The modes in this 
region have somewhat enhanced amplitudes in the region 
just outside the convective core. In PULSE and LOSC this 
apparently happens for modes somewhat different from the 



remaining codes. This is the origin of the frequency differ- 
ences. As the mode trapping in this region for these stel- 
lar models is related to the Brunt- Vaisala frequency in the 
jj. -gradient zone, the previously mentioned inadequate treat- 
ment of this region in model M2k is likely the reason for 
these differences. 

The top panel of Fig. [9] presents inter-comparisons for 
1 = 2 frequencies calculated with a second-order scheme 
without Richardson extrapolation and for model M4k, for 
the same frequency range as in Fig. [7] POSC is again ex- 
tremely close to GraCo, within 2 x 10~ 4 jUHz. FILOU and 
ADIPLS present very similar results, with NOSC show- 
ing slightly larger differences, ranging from 0.0005 up to 
0.0015/iHz relative to this group. The differences decrease 
globally as far as the frequency decreases until reaching a 
magnitude of 0.0015 fi Hz for the smallest frequency studied 
here. 

The middle panel of Fig. [9] shows the same comparison 
for codes using a fourth-order scheme or a second-order plus 
Richardson extrapolation. Once again it can be seen that the 
precision resembles the results in the previous panel. Here 
POSC and OSCROX provide very similar results to GraCo. 
ADIPLS-LOSC-NOSC present small differences among them, 
with NOSC having a small increasing difference with fre- 
quency with respect to the other two. Again the general pat- 
tern here, and in the top panel, largely reflects the differences 
in G. In this case we can distinguish the codes using fourth- 
order integration scheme with its apparently noisy profiles 
as compared with a solution using second-order integration 
plus Richardson extrapolation. 

In the bottom panel of Fig. [9] an inter-comparison for 
the different "degrees of freedom" using only GraCo is pre- 
sented; as in earlier corresponding plots the same value of 
G is used as in the computation of the equilibrium model. 
As expected the Richardson extrapolation has a growing in- 
fluence as the frequency decreases and the mode order in- 
creases, with quite substantial differences, compared with 
the differences between different codes, for the lowest- 
frequency modes. Thus, with the mesh provided by the evo- 
lution calculation the second-order schemes have inadequate 
numerical precision. In this case, the differences provided by 
the use of r/P as integration variable are negligible; The use 
of the Lagrangian perturbation to the pressure ( 8P) gives rise 
to frequency differences exceeding those obtained between 
the different codes, at the lowest frequencies; we note, how- 
ever, that a corresponding comparison between ADIPLS and 
LOSC does not show this effect which may therefore be par- 
ticular to the GraCo implementation. 

Finally, we want to point out that in the case of model 
M4k, the large differences appearing in the mode trapping 
region are not found. Fig. [10] presents the same period sep- 
aration as Fig. [8]but for this model. All the codes give quite 
similar results. Two mode trapping regions appear with the 
same frequency domain as in Fig. [8] As this model does 
not present any numerical imprecision in the Brunt- Vasala 
frequency pattern, the obvious conclusion is that, as for the 
mixed modes in avoided crossing, any numerical impreci- 
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Fig. 9 Frequency comparison (reference line is GraCo) for modes with 
I = 2 as a function of the frequency obtained for the model M4k in 
the low-frequency region. In the top panel the models with: second or- 
der, no Richardson extrapolation, and r, are depicted (NOSC-ADIPLS- 
GraCo-FILOU-POSC). The middle panel presents the differences ob- 
tained for the models providing fourth-order integration solutions or 
second-order plus Richardson extrapolation (LOSC-OSCROX-NOSC- 
ADIPLS-GraCo-POSC). In the bottom panel an inter-comparison for 
the different degrees of freedom only using GraCo is presented. 
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Fig. 10 Period spacing as a function of the frequency in the asymptotic 
g-mode region for the model M4k. 



sion in the description of N 2 coming from the equilibrium 
model can give rise to large differences in the frequencies 
calculated by different oscillation codes for the g-modes 
trapped in the fi -gradient zone. 



5 Large separations (LS). 

This section is devoted to the asymptotic behaviour of p- 
modes through the use of the so called "large separations", 
that is, the difference between two consecutive modes with 
the same spherical degree £ (A = v(n,£) — v(n — l,£), n be- 
ing the radial order of the mode). The structure of the section 
is similar to the previous one. We will use the same defini- 
tions as in the previous section to study all the frequencies 
ranges. From now on, the results with the M2k model will 
not be discussed, since no additional information is found 
from the further inter-comparisons. 



5.1 Large separation of radial modes. 

Figure QT| shows the results obtained for model M4k. In the 
top panel differences for the different codes using second 
order, no Richardson extrapolation, r, are depicted. POSC 
and GraCo are again extremely close. The rest of the codes 
(ADIPLS-NOSC-FILOU-GraCo) present differences of up to 
around 0.0015 /iHz, generally sharing an oscillatory pattern, 
particularly at relatively low frequency; we have no expla- 
nation for this behaviour. However, the effect is evidently 
small. 

In the middle panel of Fig.Qj]the LS differences for the 
different codes, using a fourth-order integration scheme or a 
second-order plus Richardson extrapolation, are presented. 
With the exception of LOSC the global behaviour is similar to 
that obtained without Richardson extrapolation, the NOSC- 
ADIPLS-OSCROX-GraCo differences being always lower than 
0.002 jxliz. The pattern of the LOSC differences, presenting 
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Fig. 11 Large separation inter-comparison (reference line is GraCo) 
for modes with I = as a function of the frequency calculated for the 
model M4k. In the top panel the models with second order, no Richard- 
son extrapolation, r, are depicted (NOSC-ADIPLS-GraCo-FILOU- 
POSC). The middle panel presents the differences obtained for the 
models using fourth-order integration solutions or second-order plus 
Richardson extrapolation (LOSC-OSCROX-NOSC-ADIPLS-GraCo- 
POSC). In the bottom panel an inter-comparison for the different "de- 
grees of freedom" using only GraCo is presented. 



differences one order of magnitude larger than the rest of 
the codes, is clearly related to the corresponding oscillatory 
pattern found in Fig. |H as before it is explained by the use 
of the LAWE differential equation. 

The bottom panel of Fig. QT| shows the LS differences for 
the radial modes obtained with GraCo and model M4k when 
different "degrees of freedom" are used. The Richardson ex- 
trapolation introduces differences increasing with frequency 
and mode order, as expected, giving the largest differences, 
around 0.07 fiHz. This again emphasizes the inadequacy of 
the second-order schemes for the highest-order modes, on 
the M4k mesh. The integration variable r/P gives differ- 
ences slightly lower than 5 x 10 -4 /iHz, i.e., much smaller 
than that found for different codes using the same numerical 
techniques, as depicted in the previous panels. The use of the 
Lagrangian perturbation to the pressure 8P and the LAWE 
differential equation show the same oscillating behaviour 
and values as those previously observed for the LOSC re- 
sults. As discussed above, this is related to the inconsistency 
in A* in model M4k (see Sect.[ 



5.2 Non-radial modes with 1 = 2 

To illustrate the differences appearing in the case of non- 
radial modes, the spherical degree 1 = 2 has been arbitrarily 
chosen. We have divided the frequency spectrum in three re- 
gions, like in the direct frequency inter-comparison: 1) high- 
frequency region, 2) frequencies around the fundamental ra- 
dial mode and 3) low-frequency region. In all cases, a study 
similar to that developed in the radial case has been carried 
out. In the low-frequency region, the more physical period 
separation is studied, instead of the frequency separation rel- 
evant for acoustic modes. 

In the first region the results are very similar to those 
obtained for the radial case; therefore the plots are not pre- 
sented here. As in the direct frequency inter-comparison 
case, LOSC also presents an oscillating pattern, owing to the 
use of 8P as eigenfunction (see Sect. 17.21) . On the other 
hand, the only noticeable difference, when compared with 
the radial case in this region, is that for 1 = 2 the preci- 
sion among codes using the same integration procedures is 
slightly higher. 

The results obtained with the codes using a second-order 
scheme, for the modes around the fundamental radial mode, 
are depicted in the top panel of Fig. [T2l POSC remains very 
close to GraCo, with the other set of codes (F1LOU-NOSC- 
ADIPLS) extending the oscillatory pattern in the top panel of 
Fig. QT1 with the largest difference being 0.01 /iHz for the 
mixed modes. This set of codes agrees to within differences 
around 0.001 jliHz (slightly larger for NOSC). 

The middle panel of Fig. [12] depicts the differences 
for a fourth-order integration scheme or second-order plus 
Richardson extrapolation. The precision of the different 
codes is similar to that given by the second-order scheme 
without Richardson extrapolation. The maximum difference 
is lower than 0.08 /iHz, and most of the codes present differ- 
ences around 0.005 jiHz. Once again the largest differences 
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Fig. 12 Large separation inter-comparison (reference line is GraCo) 
for modes with I = 2 as a function of the frequency calculated for the 
model M4k and for modes around the fundamental radial one. In the 
top panel the models with second order, no Richardson extrapolation, 
r, are depicted (NOSC-ADIPLS-GraCo-FILOU-POSC). The middle 
panel presents the differences obtained for the models using fourth- 
order integration solutions or second-order plus Richardson extrapola- 
tion (LOSC-OSCROX-NOSC-ADIPLS-GraCo-POSC). In the bottom 
panel an inter-comparison for different "degrees of freedom" using 
only GraCo is presented. 



are obtained in the mixed modes. No different behaviours 
depending on the integration scheme are found. 

The bottom panel of Fig.[l2]is devoted to the differences 
obtained with GraCo when different options for the solution 
of the differential equations are chosen. In this frequency 
region the replacement of r by r/P as integration variable 
introduces the smallest differences, one order of magnitude 
lower than those obtained using the codes with the same "de- 
gree of freedom". Using the Lagrangian perturbation to the 
pressure (8P) as eigenfunction causes some difference indi- 
cating sensitivity to whether or not A* is used. The Richard- 
son extrapolation introduces differences lower than those 
found for the different codes, except in the avoided cross- 
ing zone of mixed modes, where changes comparable to the 
largest one are found. Given the rather substantial variations 
in the bottom panel, and the fact that most codes show the 
same variation in the top two panels, one might suspect that 
the dominant source of this variation is in fact in the GraCo 
results used as reference. 

The low-frequency region is studied through the period 
separation for I = 2 (ATI in seconds), illustrated in Fig. 1131 
The top panel of this figure shows the differences found 
using only codes with a second-order integration scheme. 
ADIPLS and NOSC show similar shifts of around 0.1 s rela- 
tive to GraCo. FILOU is similar at the higher frequencies but 
shows a small oscillating behaviour in the mode-trapping re- 
gions. Finally POSC presents a quite noisy pattern, varying 
around zero. Again the overall grouping of the differences 
(ADIPLS-NOSC-FILOU and GraCo-POSC) reflects the different 
values of G. 

Results obtained using a fourth-order integration scheme 
or a second-order plus Richardson extrapolation are com- 
pared in the middle panel of Fig. Q2J The values of the dif- 
ferences found in this case are of the same order as in the 
previous inter-comparison. Here we can distinguish codes 
using a fourth-order scheme or a second-order plus Richard- 
son extrapolation, because of the apparently random pat- 
tern they present, with differences one order of magnitude 
larger than the main values. The frequency differences ob- 
tained with the OSCROX and LOSC results are those present- 
ing a noisy behaviour, when comparing with a second-order 
plus Richardson extrapolation solution (as GraCo does) as the 
reference line. POSC also presents some differences in the 
mode-trapped regions when compared with other codes but 
using the same integration scheme. 

The bottom panel of Fig. [T5] shows the differences ob- 
tained with the same code (GraCo) and different choices of 
the "degrees of freedom". The use of the Richardson ex- 
trapolation introduces substantial differences, of the order of 
seconds, with noticeable wiggles in the two mode-trapping 
regions. As noted previously this reflects the inadequacy of 
the second-order schemes for high-order modes. The use of 
the Lagrangian perturbation to the pressure as variable (dP) 
also introduces substantial differences, particularly around 
the trapped modes near 22 /iHz, related to the frequency dif- 
ferences found with GraCo in this region when 8P is used 
(cf. Fig. |9). Using r/P as integration variable gives a small 
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Fig. 13 Period separation inter-comparison (reference line is GraCo) 
for modes with I = 2 as a function of the frequency calculated for the 
model M4k in the low-frequency region. In the top panel the mod- 
els with second order, no Richardson extrapolation, r, are depicted 
(NOSC-ADIPLS-GraCo-FILOU-POSC). The middle panel presents 
the differences obtained for models using fourth-order integration 
solutions or second-order plus Richardson extrapolation (LOSC- 
OSCROX-NOSC-ADIPLS-GraCo-POSC). In the bottom panel an 
inter-comparison for the different choices of "degrees of freedom" us- 
ing only GraCo is presented. 



difference when compared with r, one order of magnitude 
lower than the differences found between the codes using 
the same numerical integration schemes. 



6 Small separations (SS). 

In this section we study the inter-comparisons for the small 
separations (SS) 8(f = v(n,£) - v(n—l,£+2). Therefore, 
two sets of inter-comparisons can be done, one for the 
(£ = 0-1 = 2) modes and another for the (£ = 1 - £ = 3) 
modes. In this case only the high-frequency region is stud- 
ied, since this is where this quantity has physical meaning. 
As in the previous section we concentrate on results for 
model M4k. 



6. 1 Small separations 5q2 

The results obtained for the small separation 8q2 are pre- 
sented in Fig. Q41 The top panel of this figure shows the 
differences obtained with the codes solving the set of dif- 
ferential equations with a second-order integration scheme. 
In the avoided crossing region some wiggles occur related 
to variations for 1 = 2, but these wiggles are of the order of 
magnitude of the differences obtained for the high-frequency 
region. ADIPLS-POSC-GraCo present similar values for high 
frequencies and NOSC-FILOU presents differences around 
0.004 jJ.Hz, and always lower than 0.01 juHz. These differ- 
ences are much lower than the expected observational errors 
for C0R0T. 

The middle panel shows the same inter-comparison us- 
ing a fourth-order integration scheme or a second-order plus 
Richardson extrapolation. Once again we cannot distinguish 
the integration scheme used. Just LOSC shows an oscillatory 
pattern for high frequencies due to the use of LAWE or dP 
(see Sect. 17.21) . The order of magnitude of the main differ- 
ences is the same as those obtained using only a second- 
order scheme. The wiggles in the avoided crossing regions 
are still present, and for high frequencies the differences 
are all in the range [-0.001,0.001] jtiHz. ADIPLS and POSC 
present an almost constant difference with GraCo, NOSC and 
OSCROX show a small noisy behaviour. 

The bottom panel of Fig. Q4] presents the SS differences 
induced in GraCo when different choices of "degrees of free- 
dom" are selected. In this case they are all below, or of the 
same order of magnitude as, the spread between the different 
codes illustrated in the middle panel. Interestingly, Richard- 
son extrapolation introduces differences for higher frequen- 
cies far smaller than found for the large separation (Fig. II lb. 
The use of the Lagrangian perturbation to the pressure (8P) 
as variable results in the same oscillatory pattern as seen for 
LOSC in the middle panel. The integration variable r/P gives 
increasing differences in the range [0.001,0.003] /iHz, sim- 
ilar to the Richardson extrapolation, probably reflecting a 
difference in the sensitivity of radial and nonradial modes to 
the choice of independent variable in the GraCo code; how- 
ever, the effect is obviously small. Note that in this case the 
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Fig. 14 Small separation inter-comparison (reference line is GraCo) 
for modes with i = — 2 as a function of the frequency cal- 
culated for the model M4k in the high-frequency region. In the 
top panel models with second order, no Richardson extrapolation, 
r, are depicted (NOSC-ADIPLS-GraCo-FILOU-POSC). The mid- 
dle panel presents the differences obtained for models using fourth- 
order integration schemes or second-order plus Richardson extrapola- 
tion (LOSC-OSCROX-NOSC-ADIPLS-GraCo-POSC). In the bottom 
panel an inter-comparison for different "degrees of freedom" only us- 
ing GraCo is presented. 



use of the LAWE for the radial modes, keeping the default 
"degrees of freedom" (including the use of P') for i = 2, is 
hardly meaningful; thus no results are included for LAWE. 

6.2 Small separations <Si3 

In Fig. Q3] the inter-comparison of 5i3 is done for model 
M4k. The top panel shows the differences given by the 
codes using a second-order integration scheme. For POSC- 
ADIPLS-FILOU-GraCo the differences are between -0.001 
and 0.002 fiHz and NOSC shows somewhat larger values. 
In any case the precision is good and the patterns rather 
smooth. 

The codes using a fourth-order integration scheme or a 
second-order plus Richardson extrapolation have produced 
the results presented in the middle panel of this figure. The 
precision is very similar as in the previous panel or even 
a little higher. LOSC shows the oscillatory pattern already 
obtained in all the previous inter-comparisons. POSC and 
ADIPLS present smooth difference profiles when compared 
with GraCo. On the other hand OSCROX and NOSC give a low 
noisy profile. Once again the integration scheme used cannot 
be discriminated. 

Finally, the bottom panel presents the effect in 5b 
of using different numerical integration schemes. Using 
Richardson extrapolation gives large differences in the high- 
frequency region, even larger than the precision of the differ- 
ent codes while, as expected, the effect is small for low-order 
modes. If the differential equations are solved with r/P as 
integration variable, an almost increased difference is intro- 
duced, always lower than 0.001 fx Hz, i.e., of the order of the 
precision among most of the codes using the same numeri- 
cal techniques. The use of the Lagrangian perturbation to the 
pressure as variable introduces the same oscillatory pattern 
as was obtained with LOSC. 



7 Computational variations 

7.1 The influence of the gravitational constant G 

As indicated in Table|2]different values of G are used by dif- 
ferent codes. In many cases the value differs from the value, 
G = 6.67232 x 10~ 8 cgs, which was used in the computation 
of model M4k. As a result, as seen by the pulsation code the 
equilibrium model is not strictly in hydrostatic equilibrium, 
thus potentially causing errors in the computed frequencies. 

In this section we examine the consequence of such in- 
consistencies. Thus we test the influence of the choice of 
a particular gravitational constant G, within always the re- 
cent experimental values found in the literature. The value 
of this constant is not very accurately known, and the differ- 
ent values we can find in the literature can have an impact 
on the frequency calculation of the same order or even larger 
than the differences studied here. To illustrate this, two ex- 
treme values of G found in the literature have been cho- 
sen: G, = 6.6716823 x 10~ 8 cgs (as fixed for Task 1), and 
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Fig. 15 SS inter-comparison (reference line is GraCo) for modes 
with £ = 1 — 3 as a function of the frequency obtained for model 
M4k in the high-frequency region. In the top panel the models with 
second order, no Richardson extrapolation, r, are depicted (NOSC- 
ADIPLS-GraCo-FILOU-POSC). The middle panel presents the differ- 
ences calculated for models using fourth-order integration solutions or 
second-order plus Richardson extrapolation (LOSC-OSCROX-NOSC- 
ADIPLS-GraCo-POSC). In the bottom panel an inter-comparison of 
the different "degrees of freedom" using only GraCo is presented. 
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Fig. 16 Frequency comparison of modes with £ = and 2 obtained 
with GraCo,(//, SP = 0, no Richardson extrapolation, r) and the model 
M4k when two different values of the gravitational constant G are used. 



G 2 = 6.693 x 10" 8 cgs dFixler et all l2007h . the most recent 
one, although with a quoted random error of ±0.027 x 10~ 8 
and a systematic error of ±0.021 x 10~ 8 cgs it is con- 
sistent with the previous value. The present recommended 
value of the Committee on Data for Science and Tech- 
nology (CODATA) can be found in the World Wide Web 
at physics.nist.gov/constants, and it is closer to that 
fixed in Task 1. The comparison has been carried out with 
GraCo at fixed "degrees of freedom" (P' , 8P = 0, no Richard- 
son extrapolation, r). The differences obtained using both 
values of G are shown in Fig. [16] for modes with i = and 
2. The equilibrium model used is M4k. 

Surprisingly, the differences obtained in the region 
around the fundamental radial mode are much larger than 
those obtained for the direct frequency inter-comparison for 
the different codes or for the same code using different in- 
tegration numerical schemes. As expected for radial modes 
the largest difference is obtained for the fundamental radial 
mode. The value of this difference (more than 0.3 jliHz in 
magnitude) is really considerable. In the asymptotic regions 
the influence of the value of this parameter decreases until 
reaching values lower than 0.1 fiHz. Nevertheless, these dif- 
ferences are larger than those obtained by using or not the 
Richardson extrapolation. The presence of wiggles for £ = 2 
for the mixed modes is also remarkable. 

The relative difference between G\ and Gj is around 
3 x 10 -3 . In contrast, the relative difference between, say, 
the value chosen for Task 1 and the ADIPLS value is only 
5 x 10~ 5 . Thus we might expect this difference to cause a 
frequency difference of order 5 x 10~ 3 /iHz. In fact, both 
the magnitude and shape of the differences between ADIPLS 
and GraCo in the middle panel of Fig. [4] are consistent with 
this estimate. In general, as already discussed, it appears that 
much of the differences between codes using the same nu- 
merical scheme results simply from the different values of 
G. Thus, it is obviously important to ensure that consistent 
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values of G are used in the evolution and oscillation calcula- 
tions. 



7.2 The choice of dependent variables and equations 

The most dramatic difference between the different codes 
is the oscillatory variation shown by the computations us- 
ing the LAWE or 8P, relative to the reference GraCo re- 
sults (e.g., Fig. lU. The oscillatory nature and the 'period' 
of the variation indicated that it reflected a sharp differ- 
ence in some aspect of the model, located approximately 
in the region of the second helium ionization zone. Af- 
ter various failed attempts to identify the cause of the 
variation it tr anspired that it resulted f r om a, previously 
known (e.g., iBoothrovd and Sackmannl 120031) . inconsis- 
tency in the OPAL equati on-of-state tables (the original ta- 
bles of Ro gers et al.Lfl996h used in the computation of model 
M4k with the ASTEC code. 

To understand this we consider the computation of 



1 dlnP dlnp 



Fi dlnr dlnr 



(1) 



From the point of view of stellar evolution calculation this is 
somewhat inconvenient, since the gradient of p does not ap- 
pear directly in the equations of stellar structure. Thus eval- 
uation of A* would involve numerical differentiation of p. 
A potentially more convenient formulation is obtained by 
rewriting the expression in terms of the temperature gradi- 
ent which is known from the equation of energy transport, 
given that p is a function, known from the equation of state, 
of pressure, temperature and composition which we charac- 
terize solely by the hydrogen abundance X. Thus we obtain 



dlnp 
dlnT 



(V-V, 



ad / 



P.X 



dlnP 
dlnr 



( d\np\ dlnZ 
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RT 



din; 
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here, as usual, V = dlnT/dlnP and V a( j is its adiabatic 
value. This expression still requires numerical differentia- 
tion of but in much of the model X varies little or not 
at all and hence the term in dlnX/dlnr makes a limited 
contribution. Furthermore, in the bulk of convective zones, 
with homogeneous composition, Eq.|2]accurately reflects the 
small value of A* resulting from the small value of V — V a( j 
determined, e.g., from the mixing-length treatment. 

The transformation in Eq.|2] obviously assumes that the 
thermodyna mical quantities used are mutu ally consistent. 
As noted by Boothro yd and Sackmannl d2003h . this is not the 
case o f the original OPAL tables provided by iRogers et al.l 
( 1996). The effect of this on A* is illustrated in Fig. [171 for 
model M4k, in terms of the difference between the value 

2 except if diffusion is taken into account; in that case the gradient 
of X appears directly in the equations ( see, fo r example, the description 
of ASTEC by Christensen-Dalsgaard ( 2007)) 




Fig. 17 Differences between A* as evaluated from Eq. [T] and evalu- 
ated from Eq.|2]as done in ASTEC, for model M4k, against fractional 
radius. The differences have been normalized by V = |dlnP/dlnr|. 



of A* calculated directly from Eq. [TJ using numerical dif- 
ferentiation to evaluate dlnp/dlnr|3and the value, A\ STEC , 
evaluated from Eq. [2] There are obvious systematic differ- 
encesQ concentrated in the ionization zones of hydrogen 
and helium. These would indeed affect the oscillations as 
a 'sharp feature' and hence cause an oscillatory signature in 
the frequenc ies, with the ' period' seen in Fig. B and else- 
where (e.g.. iGoughl 1 19901 : [Monteiro and Thompson! 120051 : 
iHoudek an d GoughLl2007h. 

To test further this interpretation of the results we 
have recomputed frequencies for model M4k, but replac- 
ing A^ STEC with A* computed from Eq.Q] For this modified 
model the differences between ADIPLS and LOSC are smaller 
by an almost an order of magnitude than the differences il- 
lustrated in Fig.|4]and with none of the systematic oscillatory 
character. 

We finally note that transforming between the equations 
using P' and 8P also depends on the equation of hydrostatic 
equilibrium and hence on a consistent choice of G in the os- 
cillation equations (cf. Sect. l7.lb . However, within the actual 
range of G-values used by the different codes this effect is 
smaller than the effect of the inconsistency in A*. 



7.3 The choice of independent variable 

We have found a fairly significant dependence on the G 
value used by the oscillation codes for the differences of 
using the integration variables r or r/P. Fig. 1181 shows this 
dependence in the case of I = 0. In this figure, the differ- 
ences between the use of r or r/P with the three different 
G values presented in Table 2, are depicted. All of them 
have been obtained using GraCo and M4K. The differences 

3 at this level of precision it does not matter whether dlnP/dlnr is 
calculated from the equation of hydrostatic support or through numer- 
ical differentiation. 

4 the small-scale rapid variation is probably associated with the in- 
terpolation procedure in the equation-of-state tables 
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Fig. 18 Differences between the use of r or r/P as integration variable 
as a function of the mode frequency, for radial modes. Different G 
values are used in the GraCo oscillation code: G (M4K) is the G value 
used to obtain the M4K model, G (M2K) is the G value used to obtain 
the M2K model, and G NOSC is the G value used by the NOSC code. 
M4K has been used as equilibrium model for this comparison. 

increase with frequency, and their value depend on the dif- 
ference G eq — G osc , where G x is the G value used to obtain 
the equilibrium model or the oscillation frequencies, respec- 
tively. In fact, the transformation between using r and r/P as 
independent variables uses the equation of hydrostatic equi- 
librium and hence the result is obviously sensitive to whether 
or not a consistent value of G is used. By far the smallest 
differences are obtained when the same G are used in the 
equilibrium and oscillation codes. 

Even when the same value of G is used, we have found a 
slight effect of using r/P as independe nt variab l e, as i mple- 
mented in a few of the codes (see also iProvostl (120071) ). in- 
stead of r; in Fig.[l8]the maximum difference is 0.008 /iHz, 
as also found in Fig. [4] The origin of choosing r/P seems 
to be a concern that near the surface, where r varies little 
between adjacent mesh points, rounding errors in the eval- 
uation of differences in r, when representing the differen- 
tial equations on finite-difference form, could have a signif- 
icant effect on the results. Given the rapid variation of P in 
the near-surface layers this problem is obviously avoided if 
r/P is used instead. While there might well be such prob- 
lems with using r when variables are represented in single- 
precision form (with four-byte reals) it seems unlikely to 
be a problem when using double -precision variables^ On 
the other hand, the present test has been carried out with 
models transferred in the FG0NG format, where variables are 
given with ten significant digits. Here one cannot exclude 
that rounding errors might be significant. 

To test this JC-D computed frequencies with ADIPLS us- 
ing both the original binary (double-precision) version of 
M4k and the version resulting from conversion to FG0NG 
format. The frequency differences between these two cases 
were always less than 10~ 5 /iHz, strongly indicating that the 

5 In contrast, the problem would be severe if the interior mass were 
used as independent variable. 



use of the FG0NG format is not a concern with the present 
level of accuracy. 

It is obvious that the choice of independent variable also 
affects the truncation error, i.e., the error in the representa- 
tion of the differential equations on finite-difference form. 
This can presumably be the reason for the differences found 
with GraCo when comparing the use of r and r/P (e.g., 
Fig. [T41 , although the effects are small. Further analysis is 
required to decide which is the more accurate representa- 
tion. 



8 Conclusions. 

A complete inter-comparison of the frequencies calculated 
by using a set of oscillation codes has been presented. Nine 
oscillation codes have been used: ADIPLS, FILOU, NOSC, 
LOSC, GraCo, OSCROX, PULSE, POSC and LNAWENR. The 
only free parameter for each code has been the integration 
numerical scheme. In order to ensure that the same physics 
has been imposed to all the codes two equilibrium models 
of a \.5M Q star have been supplied, with 4042 (M4k) and 
2172 (M2k) mesh points in their grids, respectively. No re- 
meshing has been allowed. The model M2k presents inad- 
equate numerical resolution in the Brunt-Vaisala frequency 
close to the boundary of the convective core. The model M4k 
does not suffer from these problems. In the present paper we 
have used these models only to show the sensitivities of the 
oscillation codes to possible numerical inaccuracies in the 
equilibrium models. A complete study of the accuracies of 
the different evolutionary codes is provided by Lebreton et 
al. in this volume. 

Several inter-comparisons have been performed: 1) Di- 
rect frequency inter-comparison, 2) large separation (LS), 3) 
small separation (SS), and 4) different experimental values 
of the gravitational constant G. 

The main conclusions can be summarized as: 

- When codes using the same numerical integration sche- 
mes are compared, the general precision obtained for 
model M4k (0.02 /iHz) is higher than the expected preci- 
sion of the observational data obtained with the CoRoT 
mission. This precision is always one order of magni- 
tude better or even more when compared with that ob- 
tained for model M2k (0.5 /iHz). The largest differences 
are obtained for the high-frequency region. 

- The use of a second-order integration scheme plus 
Richardson extrapolation or a fourth-order integration 
scheme does not improve the agreement between the 
codes when compared with the use of a second-order 
scheme. In the cases here presented it was not possible 
to distinguish between the use of a fourth-order scheme 
or Richardson extrapolation, which always give very ac- 
ceptable consistency for model M4k. 

- However, the use of a fourth-order integration scheme 
(or a second-order plus Richardson extrapolation) intro- 
duces differences larger than 0. 1 /iHz (1.5 /iHz for M2k 
and 0.8 /iHz for M4k) compared with a second-order 
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scheme, mainly located in the regions of high order (high 
frequency for p-modes and low frequency for g-modes). 
This indicates that the second-order schemes have in- 
adequate numerical precision on the mesh provided in 
M4k and, obviously, even more so on the sparser mesh 
in M2k. 

- For "large separations" the agreement between different 
comparable codes using model M4k, with every possible 
numerical scheme is generally better than the expected 
precision of the data. However, the equation-of-state in- 
consistency in model M4k (see Sect. 17.21) introduces a 
signature that it potentially significant. The use of the 
Richardson extrapolation together with a second-order 
integration scheme, or a fourth-order scheme, gives dif- 
ferences for high-order p- and g-modes relative to the 
use of only a second-order scheme that are comparable 
to the expected observational accuracies, indicating that 
the accuracy of the latter scheme is inadequate. 

- The "small separations" present differences for model 
M4k lower than the expected observational accuracies 
everywhere, although the inconsistency in M4k has a no- 
ticeable effect. Use of Richardson extrapolation has a po- 
tentially significant effect, indicating that the precision 
of second-order schemes may be inadequate. 

- The use of r/P as integration variable instead of r has 
a small influence when a consistent value of G is used. 
However, an inconsistency between the G values used to 
obtain the equilibrium model and the oscillation frequen- 
cies can give substantial differences between the use of 
these integration variables, over the range of G values in 
Table 2. 

- For model M2k, the mixed modes near the avoided 
crossing present large differences, reaching even up to 
4/iHz. These differences are not present when model 
M4k is used. The problems with TV 2 in M2k are un- 
doubtedly the main source of this difference. The same 
behaviour is found for the trapped g-modes. Therefore, 
the correct numerical description of TV 2 is critical for the 
value of the frequencies of these trapped modes or in 
avoided crossing. 

- The use of dP or P' as eigenfunction, or the solution of 
the LAWE, may have a significant effect if the equilib- 
rium model is based on thermodynamic quantities that 
are not internally consistent, as is the case for the OPAL 
tables used to calculate model M4k. Further effects arise 
when there are problems with the Brunt- Vaisala fre- 
quency (as for model M2k); here a different choice can 
change the frequencies of some modes since N 2 does 
not appear in the differential equations when 8P or the 
LAWE are used and it does when using P'. 

- The value of the gravitational constant G in the oscil- 
lation calculations can introduce non-negligible differ- 
ences as well, if it is not the same as the value used 
in the equilibrium model. When the two extreme values 
found in the literature are used, such inconsistency yields 
differences in the range [—0.35, —0.08] /iHz for model 
M4k, larger than those obtained when different numer- 



ical integration schemes are used. Differences between 
the values of G actually used by the different codes, al- 
though less extreme, account for much of the difference 
between the computed frequencies. 

Therefore, for a proper pulsational study, we require that 
the number of mesh points and their distribution must be 
such as to yield an equilibrium model that satisfies the dy- 
namical equations with sufficient accuracy also in the re- 
gions of the star where the physical quantities present rapid 
variations (e.g., the outer layers and /i-gradient zones). In 
addition, the mesh used in the pulsation calculation must 
properly resolve the eigenfunctions of the highest-order 
modes considered. In the present case model M4k, with 
4042 points, appears to satisfy these conditions although the 
use of a fourth-order integration scheme, or a second-order 
scheme and Richardson extrapolation, is still needed in the 
oscillation calculation; these higher-order schemes give sig- 
nificant improvements compared with the use of a simple 
second-order scheme. The use of a second-order integration 
plus Richardson extrapolation scheme is not distinguishable, 
in accuracy and precision, from the use of a fourth-order in- 
tegration scheme. A correct physical and numerical descrip- 
tion of the Brunt-Vaisala frequency is essential when P' is 
used as eigenfunction; in particular, inconsistencies in the 
equation of state can have serious effects on the frequencies. 
Inconsistency between the value of G used in the oscillation 
calculation and the value used to compute the equilibrium 
model, within the range of the different values of G found in 
the literature, may lead to substantial errors in the computed 
frequencies. 

We note that the situation is somewhat different if con- 
sistent values of G are used in the evolution and oscilla- 
tion calculations. Then the effect on the frequencies is ap- 
proximately given, according to homology arguments, as a 
scaling by (GM) 1 / 2 . However, since the product GM is 
known extremely precisely from planetary motion in the 
solar system, any variation in G should be reflected in a 
corresponding change in the assumed value of Mq. If this 
is the case, and if the model is characterized by a given 
value of M/Mq (as is typically the case) the effect on 
the frequencies of changes in G are very small (see also 
IChristensen-Dalsgaard et all 120051) . 

In further tests more care is required to secure the full 
consistency of the models: a consistent equation of state 
should be used, and the value of G should obviously be the 
same in the equilibrium-model and the pulsation calcula- 
tions; indeed, this strongly argues for including the value of 
G as one of the parameters in the model file. The main con- 
clusion of this extensive investigation, however, it positive: 
with a properly resolved equilibrium model the broad range 
of oscillation codes likely to be involved in the asteroseis- 
mic analysis of data from CoRoT and other major upcom- 
ing projects generally give consistent results, well within the 
expected errors of the observations. Thus, although the re- 
maining problems in the calculation evidently require atten- 
tion, we can be reasonably confident in our ability to com- 
pute frequencies of given models and hence in the inferences 
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concerning stellar structure drawn from comparing the com- 
puted frequencies with the observations. 
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